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Abstract 

We compute by Monte Carlo numerical simulations the critical exponents 
of two-dimensional scalar field theories at the tricritical point. The re- 
sults are in agreement with the Zamolodchikov conjecture based on conformal 
invariance. 
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In the last decade there has been a considerable progress in the classification of uni- 
versality classes of two-dimensional (2D) critical phenomena in terms of conformal field 
theories [1]. For most of the universality classes there are integrable models whose critical 
exponents can be exactly calculated. However, the critical behavior of other 2D models has 
to be investigated by other means because they cannot be exactly solved. Self-interacting 
scalar theories with even polynomial interaction 

-r 2 N 

^W = 2E(5»-<^)' + E^2.<^^'- (1) 

i=l r = l 

belong to this second class of models. 

The phase structure of these theories is very complicated. It contains in general l- 
critical manifolds (/ = 2, . . . , iV) of dimensionality c/; = iV-|-l — / in the space of the coupling 
constants {\2rYrZ1 ■ A.B. Zamolodchikov [2] conjectured that the leading multicritical 
behavior of these models for arbitrary N > 2 corresponds to the conformal models of the 
unitary series with central charge 

_ 6 
'"^"(iV + l)(iV + 2)- 

However, there is no direct proof of this conjecture in terms of explicit calculations of 
correlation functions or finite-size scaling (FSS) analysis of thermodynamical quantities. 

The critical exponents can be obtained from the natural identification pointed out 
by Zamolodchikov between relevant local perturbations and primary fields (pp^q of the 
corresponding conformal theory. These fields possess scaling dimensions given by the Kac 
formula 

_ [(N + 2)p-{N + l)q]'-l 
4(iV + l)(iV + 2) ' 
l<p<iV; l<g<iV + l 

The order parameter (f> (i.e. the magnetization) is identified with the most relevant 
non-trivial primary field (f> ^2,2 with scaling dimension /i2,2 = 3/[4(iV -|- l)(iV -|- 2)]. The 
critical exponent rj can be read off from the decay of the two-point correlation function 



w-)m) ^ p^i^, (4) 

at large distances \x\ ^ 1. Thus, rj = 4/^2,2 for d = 2. 

The energy density operator e is given by the renormahzed composite field : and 
it is identified with the primary field e ^3 3 (resp. ^1,3) for iV > 3 (resp. N = 2). The 
corresponding scaling dimension is hs^s = 2/[(N + l)(iV + 2)] (resp. hi^3 = 1/2). The 
two-point correlator of the energy density has a similar behavior to that of Eq. (4); but 
now the power is related to the critical exponent v [3] 

This implies that l/v = 2(1 — /i3,3) (resp. l/v = 1). The remaining exponents can be 
obtained via the well-known scaling relations [4]. The result is displayed in Table I. Notice 
that the Gaussian critical exponents are recovered in the limit iV — )• 00. 

The free scalar bosonic field (N = 1) is the only model of the family which is exactly 
solvable. The critical massless regime describes a conformal field theory with central 
charge c = 1. Thus, its critical exponents are given by the case iV = 00 in Table I. 

The first non-trivial case corresponds to iV = 2. This model contains one critical line 
in the coupling constant space with the Gaussian model at one end (A4 = A2 = 0) and the 
Ising model (A4 — )• 00, A2 — > — oo,A4/A2 — > constant) at the other one [5]. It is generally 
believed that the whole critical line (except the Gaussian endpoint) belongs to the Ising 
universality class (c = 1/2). However, some authors have claimed that this could not be 
the case [6], although their results are not very conclusive and they might be interpreted 
as strong corrections to scaling, as suggested in [7]. 

The case iV = 3 is the first one where tricritical points are found. It is remarkable 
the striking difference between the conjectured value for rj = 3/20 and the prediction 
coming from standard perturbation theory rj = 1/500 [4]. This model was also considered 
by Hamber [8], who used an analysis based on an approximate renormalization-group 
transformation on the equivalent one-dimensional quantum Hamiltonian. A tricritical 
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point was found and the leading eigenvalue was equal to yn = l/v = 1.661 (resp. 1.817) 
for block size 6 = 2 (resp. h = 3). From Table I we see that Zamolodchikov's prediction is 
l/v = 1.8, very close to the latter result. However, a complete study of this model is still 
lacking. 

In this letter the validity of this conjecture is analyzed by explicit calculation of the 
corresponding critical exponents in the case iV = 3. We consider a square L X L lattice 
with unit spacing and periodic boundary conditions. To each site i of the lattice we assign 
a continuous unbounded spin (pi. The partition function Z of this model on the lattice is 
given in terms of the discretized Euclidean action S 

Z = JV<P e-^(^); Si<P) = - ^ <P.<P, + J2Vi<Pi) (6) 

<i,j> i 

where T>(f> = Hi d^Pi] and the first sum is over all the nearest-neighbor pairs of spins. The 
one-particle potential V is defined as 

-r JV=3 

= -(4 - m V + E (7) 

r = 2 

For convenience, we have identified A2 with minus the mass squared A2 = —m^. 

The phase diagram of this model can be described qualitatively by the mean field 
approximation [4]. For any Ae > there is a curve of singular points in the (m^,A4) 
subspace. For positive A4, the transition is second-order. However, when A4 becomes 
negative, there is a point where the transition changes its character to first order. This 
is the tricritical point. Actually, there is a line of tricritical points in the whole coupling 
constant space (m^,A4,A6). It is expected that all these tricritical points belong to the 
same universality class, irrespective of the value of Ae > 0. The numerical results reported 
here correspond to the choice Ae = 1. 

We have carefully analyzed the phase diagram by means of the variational method 
introduced in [9]. It turned out to be very accurate for the search of critical points in the 
subspace (m^,A4,A6 = 1). This method also allows to distinguish clearly among critical, 
tricritical and first-order points. In practice, two different critical points has been selected. 
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One is an ordinary (Ising-like) critical point and the other one is a critical point in the 
vicinity of the tricritical point. The region around those points has been explored by 
Monte Carlo (MC) simulations, which confirmed that they are actually very close to the 
real critical points. 

To simulate the model given by Eq. (6) we have used a heat-bath algorithm based on 
one proposed in ref. [5] for the case N = 2. This one has the important drawback that it 
only works for < 4. We have generalized their algorithm for arbitrary N by adding and 
subtracting a m'l^cfP' term to the potential in Eq. (7). The probability density of choosing 
the new spins can be written as in [5]: it is a product of a Gaussian distribution (with their 

substituted by — 2m^) and an exponential term h{cf)) = ex.]){m^cfP' — A4^'* — Ae^^), 
which is bounded if Ae > 0. This latter function can be normalized by adding some 
constant to the potential (7) which is irrelevant when computing expectation values. The 
parameter is free, except for the constraint 2m^ > — 4. Thus, it can be tuned to 
maximize the performance of the algorithm (see below). 

The following quantities have been measured: energy density E = jj{S), specific heat 
Cv = jj(T^{S) (where cr(-) is the standard deviation), magnetization M = jj{\J2i4>i\), 
magnetic susceptibility x = jj'^^il J2i ^d)? mean value of <f>^ (i.e. M' = jj{J2i ^Pl)) ^i^d 
the corresponding "susceptibility" x' = ]jcr^(Z]i^D- 

Given the values of A4 and Ae = 1 (selected by the variational method), for each L 
we vary until a maximum in any of the quantities Q = {Cv^ Xt x'} attained. This 
is done in practice by using the spectral density method (SDM) [10]. For each quantity 
Q we obtain the position m^(Q) of the peak and the height of such peak Q'^"'^ . First, 
some preliminary runs were made to select the region of interest in the axis. The 
longest runs were carried out at masses such that their distance to the predicted position 
of the peaks are less than the range of validity of the SDM Am^ 2/[La/x' ] [5]. In some 
occasions several masses were selected to cover the three peaks. A summary of the runs 
is contained in Table II. 

Using standard ESS analysis [11] we can extract from the MC data the critical mass 
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ml and the critical exponent associated with each quantity Q (i.e. the exponent p defined 
by Q{m^) ~ — m^|~'') divided by v. In our case, we can estimate the ratios 7/1^, alv 
and 7'/^ (this one corresponds to the susceptibihty x')- Notice that Zamolodchikov's con- 
jecture imphes that 7' = a. FinaUy, it is also worthy to consider the quantity M{m\{x)) 
which is defined as the value of the magnetization at the point where the magnetic sus- 
ceptibility peaks. This quantity is expected to go to zero as L~^^'^ for L large enough [12]. 
The ratio fi/v is related directly to 7/1^ via the scaling relation 

= 1 - 7/(2//) (8) 

Thus, the validity of this important relation could be tested numerically in a very clean 
way for this model. 

To check the quality of our results we have estimated the autocorrelation time for the 
magnetization {tm) and the energy density using the methods of ref. [13]. It turned out 
that the slowest mode was the magnetization, although we do not exclude the existence 
of another slow mode not considered here. If we fit our results to tm = AL^ , we obtain a 
dynamical critical exponent roughly equal io z 2.0 for the two MC simulations reported 
here (see Table II). This is in agreement with the conclusions of [14]. We have also found 
that the constant A is larger for the simulation near the tricritical point. Unfortunately, 
due to the large autocorrelation time of the heat -bath algorithm, we have been unable to 
obtain very good statistics. For the largest lattice we only could do ~ 10"^ tm measures 
(resp. ~ 2 • 10"^ tm measures) for the MC simulation close to the tricritical point (resp. 
Ising-like critical point). For the smaller lattices these numbers are one order of magnitude 
larger. Furthermore, the possibility of performing a reasonable simulation on larger lattices 
is beyond of our computer facilities. Finally, to save on CPU time, we have skipped a 
certain number of MC steps between measures. This number varies from 20 in the smaller 
lattices to 2 • 10"^ in the largest one (see Table II). 

First, we have analyzed a critical point on the expected Xcf)"^ regime by fixing A4 = (and 
Ae = 1). The results of the simulations are displayed in Table III. The three extrapolated 
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critical values of ml are compatible and converge to a value ml 4.67 it 0.08. The results 
for the ratios of critical exponents are the following 

- = 1.73 ±0.02 - = 0.14 ±0.02 - = ^ = 0. (9) 

U V V V 

The null results have been obtained from the failure of any power fit (e.g. the best 
fit gives ajv ~ 0.1 with a very large value of ~ 6 with one degree of freedom). 
On the other hand, the logarithmic fit is reasonably good (x^ ~ 0.5 with one degree 
of freedom). Similar numbers are obtained for the susceptibility x'. Errors in (9) have 
been estimated as the double of the standard deviation (i.e. 95% confidence level). The 
results agree within errors with Zamolodchikov's conjectures for the \<^'^ (N = 2) critical 
point (7/1^ = 1.75, fi Iv = 0.125, ajv = 7'/ v = 0). We also remark that the scaling relation 
Eq. (8) is satisfied within statistical errors. 

We now analyze the results for the second critical point, of coordinates A4 = —4.55 
and Ae = 1 (see Table IV). This value of A4 was selected by the variational method 
mentioned above. When the mass m^ is chosen suitably, it is presumably very close to the 
real tricritical point (m^^, Xi^tc^ Xe^tc = 1) • such a case, even if we are not exactly at the 
tricritical point, we can extract information about the critical exponents of the tricritical 
regime. For not too large lattice size L, the behavior near the tricritical point is expected 
to be very similar to that at the tricritical point (see for instance [15]). That is because for 
finite L any thermodynamic quantity is an analytical function of the coupling constants, 
even at the critical points. Close to the tricritical point the behavior of thermodynamic 
quantities interpolate between the tricritical regime and the infinite volume asymptotic 
regime. If we are really very close to the tricritical point the size of the lattice where the 
crossover between the two regimes appears might be very large. However, if the size of 
the system becomes large enough, the true asymptotic behavior is always recovered. The 
window where the tricritical regime is observed increases as we approach the tricritical 
point, therefore, it is very important to make the simulation very close to it. 

The three extrapolated critical values of m^^ are compatible and converge to a value 
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m^^ —9.815 it 0.003. The ratios of critical exponents obtained from the FSS analysis 

7 /3 

- = 1.81 ±0.02; - = 0.083 ±0.005 

u V 

- = — = 1.54 ±0.02 (10) 

u V 

are in fairly agreement with the values predicted according to Zamolodchikov's conjecture, 
-ilv = 1.85, = 0.077, a/;/ = -f'/iy = 1.60. 

The exponents (10) are very different of those associated to the \(f>'^ critical behavior 
(specially the ratios a/v = 7'/ v and /3/iy) and of those characteristic of a first-order phase 
transition (i.e. ^ I v = a I v = 7'/^ = 2, = [16]). However, all of them except fi/v are 
slightly below the conjectured results and the differences are several standard deviations 
and thus, they are statistically significant. This fact can be understood if our point is a 
little bit away from the true tricritical point, on the line of Ising-like critical points (i.e. 
A4 > A4^tc) and the lattices sizes are not too large. On the other hand, the value of is 
slightly higher than the predicted one. This is in agreement with the fact that the scaling 
relation (8) is also verified. In this context, it would be very interesting to study in greater 
detail the crossover from Ising-like behavior to first-order transitions passing through the 
tricritical point. 

In summary, our results confirm the validity of Zamolodchikov conjecture for the case 
iV = 3, and in other words, the power of conformal field theory methods into the descrip- 
tion of the critical behavior of 2D field theories. Our analysis makes more plausible that 
the conjecture holds for the whole series of XcfP''^ multicritical points, although a rigorous 
proof is still lacking. 
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TABLES 

TABLE L Critical exponents predicted by conformal invariance according to Zamolodchikov's 
conjecture. 



Exponents 



N 



N 



A > 3 



A = 00 



a 

7] 

S 



1 

4 

1 

8 

15 



37 
36 



_3_ 

20 



24 



77 

3 



2 ^ (JV+l)(JV+2)-2 

2(JV+l)(JV+2)-3 
2(JV+l)(JV+2)-4 

(JV+l)(JV+2)-4 
(JV+l)(JV+2)-2 

3 

(JV+l)(JV+2) 

3 

4(JV+l)(JV+2)-8 

4(JV+l)(JV+2) _ -, 
3 ^ 



1 

2 

1 
1 




00 
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TABLE II. Monte Carlo simulations performed in the coupling constant space (m^, A4, Xq = 1). 
We also display the size of the lattice and the total number of measures and iterations. 

A4 L w? ^ measures ^ iterations 

8 4.20 10^ 2 • 10^ 

8 4.70 10^ 2 • 10^ 

16 4.40 10^ 2 • 10^ 

16 4.60 10^ 2 • 10^ 

32 4.50 4-10'* 2 • 10^ 

32 4.60 4-10'* 2 • 10^ 

64 4.55 4-10'* 4 • 10^ 

64 4.60 4-10'* 4 • 10^ 

-4.55 8 -9.80 4-10'* 2-10^ 

-4.55 16 -9.81 10"* 2-10^ 

-4.55 16 -9.82 10"* 2 • 10^ 

-4.55 32 -9.816 1.05-10'* 1.05-10^ 

-4.55 64 -9.815 1.45-10'* 2.9-10^ 
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TABLE III. For the simulation performed at A4 = and Ae = 1, we show the values of the 
maximum of the specific heat and its position (C^l^ , m'ji^(Cy )) as functions of the lattice size L. The 
same values for the magnetic susceptibility (xl "^7 ''^i(x)) the x' susceptibility {x'l"'^-, ''^i(x')) 
are also displayed, as well as the magnetization M{w?j^{x)) at the pseudo-critical point w?j^{x)- The 
parentheses denote the errors (one standard deviation, ~ 67% confidence level) in the last relevant 
figures. 



L 


8 


16 


32 


64 




1.7637(97) 


2.228(13) 


2.675(31) 


3.120(30) 




4.6857(83) 


4.6079(60) 


4.5952(77) 


4.6176(37) 


Xl 


2.291(10) 


7.458(51) 


24.95(32) 


82.3(12) 




4.1844(58) 


4.3817(35) 


4.5039(36) 


4.5712(27) 


Xl 


0.3684(18) 


0.4327(25) 


0.4927(49) 


0.5462(47) 


^lix') 


4.512(13) 


4.5576(93) 


4.5863(87) 


4.6094(47) 


Mimlix)) 


0.4183(12) 


0.37205(86) 


0.3382(22) 


0.3059(22) 


TABLE IV 


. The same as 


Table III but for the values A4 = —4.55 and Xq 


= 1. 


L 


8 


16 


32 


64 




4.368(24) 


9.62(12) 


25.18(26) 


69.77(48) 




-9.80742(65) 


-9.81543(76) 


-9.81518(16) 


-9.81485(5) 


Xl 


16.603(71) 


59.78(39) 


212.9(16) 


739.2(51) 


^lix) 


-9.82953(84) 


-9.81843(76) 


-9.81577(17) 


-9.81495(6) 


Xl 


27.78(10) 


83.4(71) 


247.8(20) 


717.1(52) 


mlix') 


-9.82751(83) 


-9.81800(75) 


-9.81565(17) 


-9.81492(5) 


Mimlix)) 


0.69989(85) 


0.6555(16) 


0.6176(12) 


0.5842(15) 
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